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Abstract. We consider general integrable systems on graphs as discrete flat connections with the values 
in loop groups. We argue that a certain class of graphs is of a special importance in this respect, namely 
quad-graphs, the cellular decompositions of oriented surfaces with all two-cells being quadrilateral. We 
establish a relation between integrable systems on quad-graphs and discrete systems of the Toda type on 
graphs. We propose a simple and general procedure for deriving discrete zero curvature representations 
for integrable systems on quad-graphs, based on the principle of the three-dimensional consistency. Thus, 
finding a zero curvature representation is put on an algorithmic basis and does not rely on the guesswork 
anymore. Several examples of integrable systems on quad-graphs are considered in detail, their geometric 
£^ ' interpretation is given in terms of circle patterns. 

1 Introduction 

J> \ Discrete (lattice) systems constitute a well-established part of the theory of integrable systems. 



They came up already in the early days of the theory, see, e.g. [HI], [H2], and took gradually more 



and more important place in it, cf. a review in [NCI. Nowadays many experts in the field agree 
that discrete integrable systems are in many respects even more fundamental than the continuous 
ones. They play a prominent role in various applications of integrable systems such as discrete 



differential geometry, see, e.g., a review in [ BP3 ] 



Traditionally, independent variables of discrete integrable systems are considered as belonging 
to a regular square lattice I? (or its multi-dimensional analogs Z d ). Only very recently, there 
appeared first hints on the existence of a rich and meaningful theory of integrable systems on 
. non-square lattices and, more generally, on arbitrary graphs. The relevant publications are almost 

exhausted by |NDj , flOPf , |A|, [|NS]], [|KN2| , [[BH|], |Bl|, [|A| 



We define integrable systems on graphs as flat connections with the values in loop groups. This 
is is very natural definition, and experts in discrete integrable systems will not only immediately 
accept it, but might even consider it trivial. Nevertheless, it crystallized only very recently, and 
seems not to appear in the literature before |BHE ] , |BH|, |A3| . (It should be noted that a different 



framework for integrable systems on graphs is being developed by S.P.Novikov with collaborators 

i, m, in.) 



We were led to considering such systems by our (with T. Hoffmann) investigations of circle 
patterns as objects of discrete complex analysis: in Refs. pHSfl , pH| we demonstrated that certain 
classes of circle patterns with the combinatorics of regular hexagonal lattice are directly related to 
integrable systems. To realize this, it was necessary to formulate the notion of an integrable system 
on a graph which is not a square lattice (the regular hexagonal lattice is just one example of such 
a graph). 

Another context where such systems appear naturally is the differential geometry of discretized 
surfaces, see [BP3]. It was realized in the latter reference that natural domains where discrete 
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analogs of parametrized surfaces are defined are constituted by quad-graphs, i.e. cellular decomposi- 
tions of surfaces whose two-cells are all quadrilaterals (see Fig.l in [ BP3[| ). However, a systematical 
study of integrable systems on quad-graphs was not put forward in [BP3|. 

Finally, one of the sources of the present development is the theory of integrable discretizations 
of one-dimensional lattice systems of the relativistic Toda type Q. It was demonstrated by V. Adler 
| A2 | that these discretizations are naturally regarded as integrable systems on the regular triangular 
lattice, which eventually led him to the discovery of systems of the Toda type on arbitrary graphs 

In the present paper we put forward the the viewpoint according to which integrable systems on 
quad-graphs are very fundamental and explain many known phenomena in the field. In particular, 
we shall present the derivation of zero curvature representations with a spectral parameter for 
such systems from a very simple and fundamental principle, namely the consistency on a three- 
dimensional quad-graph. It should be noted that usually the zero curvature (or Lax) representation 
is considered as a rather transcendental attribute of a given integrable system, whose existence is 
easy to verify but very difficult to guess just by looking at the equation. Our approach shows how to 
derive a zero curvature representation for a given discrete system provided it possess the property of 
the three-dimensional consistency. Taking into account that, as widely believed, the whole universe 
of two-dimensional continuous integrable systems (described by differential-difference or partial 
differential equations) may be obtained from the discrete ones via suitable limit procedures, we 
come to the conclusion that there is nothing mysterious in zero curvature representations anymore: 
it can be found in a nearly algorithmic way. Also, we solve a lesser mystery of the existence of Toda 
type systems on arbitrary graphs |A3| , in that we trace there origin back to much more simple and 
fundamental systems on quad-graphs. 

The structure of the paper is as follows: in Sect. |2| we present the definition of integrable systems 
on graphs in more details. In Sect. |3| we consider quad-graphs and discuss their relation to more 
general graphs. Sect. |] contains generalities on integrable systems on quad- graphs (four-point 
equations) and their relation to the discrete Toda type systems (equations on stars). Sect. |5]-|7] are 
devoted to a detailed account of several concrete integrable systems on quad-graphs, namely the 
cross-ratio system and two of its generalizations. This account includes the algorithmic derivation 
of the zero curvature representation, discussion of the duality and related constructions, derivation 
of the corresponding systems of the Toda type along with their zero curvature representations, and 
a geometric interpretation in terms of circle patterns. Finally, conclusions and perspectives for the 
further research are given in Sect. H 



2 Discrete flat connections on graphs 

We consider "integrable systems on graphs" as flat connections with the values in loop groups. 
More precisely, this notion includes the following ingredients: 

• A cellular decomposition Q of an oriented surface. The set of its vertices will be denoted 
V(Q), the set of its edges will be denoted E(Q), and the set of its faces will be denoted F{Q). 
For each edge, one of its possible orientations is fixed. 

• A loop group G[X], whose elements are functions from C into some group G. The complex 
argument A of these functions is known in the theory of integrable systems as the "spectral 
parameter" . 

• A "wave function" ^ : V(Q) i— > G[\], defined on the vertices of Q. 

• A collection of "transition matrices" L : E{Q) i— > G[X] defined on the edges of Q. 
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It is supposed that for any oriented edge e = (v\,v 2 ) G F(Q) the values of the wave functions in its 
ends are connected via 



fKA) = L(e,A)%,A). (2.1) 

Therefore the following discrete zero curvature condition is supposed to be satisfied. Consider 
any closed contour consisting of a finite number of edges of Q: 

ei = (vt,v 2 ), c 2 = (v 2 ,v 3 ), e n = (v n ,vi). 

Then 

L(e n ,A)---L(e 2 ,A)L(e!,A) = /. (2.2) 
In particular, for any edge e = (vi,v 2 ), if e _1 = (v 2 ,vi), then 

L(e" 1 ,A)= (L(e,A)) _1 . (2.3) 



Actually, in applications the matrices L(c, A) depend also on a point of some set X (the "phase 
space" of an integrable system), so that some elements x(e) G X are attached to the edges e of 
Q. In this case the discrete zero curvature condition ( |2.2| ) becomes equivalent to the collection of 
equations relating the fields x(ei), . . . , x(e p ) attached to the edges of each closed contour. We say 
that this collection of equations admits a zero curvature representation. 

For an arbitrary graph, the analytical consequences of the zero curvature representation for a 
given collection of equations are not clear. However, in case of regular graphs, like those generated 
by the square lattice Z + iL C C, or by the regular triangular lattice Z + e 2m ^Z C C, such 
representation may be used to determine conserved quantities for suitably defined Cauchy problems, 
as well as to apply powerful analytical methods for finding concrete solutions. 



3 Quad-graphs 

Our point in this paper is that, although one can consider integrable systems on very different kinds 
of graphs, there is one kind - quad-graphs - supporting the most fundamental integrable systems. 

Definition 1 A cellular decomposition Q of an oriented surface is called a quad-graph, if all its 
faces are quadrilateral. 

Before we turn to illustrating this claim by concrete examples, we would like to give a construc- 
tion which produces from an arbitrary cellular decomposition a certain quad-graph. Towards this 
aim, we first recall the notion of the dual graph, or, more precisely, of the dual cellular decomposi- 
tion Q* . The vertices from V{Q*) are in one-to-one correspondence to the faces from F(Q) (actually, 
they can be chosen as some points inside the corresponding faces, cf. Fig. [J). Each e G E(Q) sepa- 
rates two faces from F(Q), which in turn correspond to two vertices from V(Q*). A path between 
these two vertices is then declared as an edge e* G E{Q*) dual to e. Finally, the faces from F{Q*) are 
in a one-to-one correspondence with the vertices from V(Q): if vq G V(G), and t»i , . . . ,v n G V(Q) 
are its neighbors connected with v$ by the edges ei = (vq, vi), . . . t n = (vo,v n ) G E(Q), then the 
face from F{Q*) corresponding to vq is defined by its boundary ej U . . . U e* (cf. Fig. |2|). 



Now, following [H, we introduce a new complex, the double T>, constructed from Q, Q* (notice 
that our terminology here is different from that of [jMj: the complex T> is called a "diamond" 
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Figure 1: The vertex from V(Q*) dual to 
the face from F{Q). 



Figure 2: The face from F(Q*) dual to the 
vertex from V{Q). 



there, while the term "double" is used for another object). The set of vertices of the double T> is 
V(V) = V{G) U V(G*). Each pair of dual edges, say e = (v u v 2 ) and e* = (/i,/ 2 ), as on Fig.g 
defines a quadrilateral (vi, fi, v%, /2), and all these quadrilaterals constitute the faces of the cell 
decomposition (quad-graph) V. Let us stress that the edges of T> belong neither to E{Q) nor to 



Quad-graphs T> coming as doubles have the following property: the set V(T>) may be decom- 
posed into two complementary halves, V(T>) = V(Q) U V{Q*) ("black" and "white" vertices), such 
that two ends of each edge from E(T>) are of different colours. Equivalently, any closed loop 
consisting of edges of T> has an even length. We shall call such quad-grahs even. 

Conversely, any even quad-graph V may be considered as a double of some cellular decompo- 
sition Q. The edges from E(Q), say, are defined then as paths joining two "black" vertices of each 
face from F(T>). (This decomposition of V(T>) into V{Q) and V{Q*) is unique, up to interchanging 
the roles of Q and Q* .) 

Notice that if T> is not even (i.e. if it admits loops consisting of an odd number of edges), 
then one can easily produce from T> a new even quad-graph V, simply by refining each of the 
quadrilaterals from F{T>) into four smaller ones. 

In the present paper we are intereted mainly in the local theory of integrable systems of quad- 
graphs. Therefore, in order to avoid global considerations, we always assume (without mentioning 
it explicitly) that our quad-graphs are cellular decompositions of a topological disc (open or closed). 



E{g*). See Fig.|. 




4 



In particular, our quad-graphs T> are always even, so that Q and Q* are well-defined. 



We shall consider in this paper several integrable systems on quad-graphs. For all these systems 
the "fields" z : V(T>) i— > C will be attached to the vertices rather than to the edges of the graph, so 
that for an edge e = (vq,vi) € E(T>) the transition matrix will be written as 

L(e, A) = L(zi,z 0l A). 

Here and below z^ = z{vk). It is easy to understand that the condition ( fO| ) for the transition 
matrices reads in the present set-up as 

£(zo,2i,A) = (l(zi,z ,X)) , (3.1) 
while the discrete zero curvature condition ( |2.2| ) in the case of quad-graphs is equivalent to 

L(z 2 , zi,X)L(zi,z , A) = L(z 2 , z 3 , X)L(z 3 , z , A). (3.2) 
In all our examples this matrix equation will be equivalent to a single scalar equation 

$(z , z 1} z 2 , z 3 ) = 0, (3.3) 

relating four fields sitting on the four vertices of an arbitrary face from F(T>). Moreover, in all our 
examples it will be possible to uniquely solve the equation (4.3) for any field Zq, . . . , z 3 in terms of 
other three ones. 

4 Equations on quad-graphs and discrete systems of the Toda type 

Of course, the above constructions might seem well-known, and there are several four-point inte- 
grable equations of the type (|i~3|) available in the literature, see, e.g., [HI] , [fig] , jNQC |, |Ny], |FV| , 



| BP1 |, [BP2], Al ] . However, these equations were always considered only on the simplest possible 
quad-graph, namely on the square lattice V(T>) = 1? . We will see, however, that the novel feature 
introduced in the present paper, namely considering quad-graphs of arbitrary combinatorics, has 
interesting and far reaching consequences. One of such consequences is the derivation of equations 



of the discrete Toda type on arbitrary graphs which were first introduced in [A3]. 

At this point we have to specify the set-up further. We shall assume that the transition matrices 
depend additionally on a parameter a = a(e) assigned to each (non-oriented) edge: 

L(e, A) = L(z 1 ,z ,a,X). 

So, actually, an integrable system will be characterized by a function a : E(D) i— > C. For some 
reasons which will become clear soon, the function a will always satisfy an additional condition. 

Definition 2 A labelling of E(D) is a function a : E(T>) i— > C satisfying the following condition: 
the values of a on two opposite edges of any quadrilateral from F(D) are equal to one another. 

This is illustrated on Fig.^. Obviously, there exist infinitely many labellings, all of them may be 
constructed as follows: choose some value of a for an arbitrary edge of T>, and assign consecutively 
the same value to all "parallel" edges along a strip of quadrilaterals, according to the definition of 
labelling. After that, take an arbitrary edge still without a label, choose some value of a for it, 
and extend the same value along the corresponding strip of quadrilaterals. Proceed similarly, till 
all edges of T> are exhausted. 
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Figure 4: A face of the labelled quad-graph 

The condition Q3.1| ) for the transition matrices reads in the present set-up as 

L(z ,zi,a, A) = (L(zi,zo,a, X)) , (4.1) 

while the discrete zero curvature condition (|3.2| ) in the notations of Fig. |] reads as 

L(z 2 ,Zi,a 2 ,\)L(zi,z ,ai,\) = L(z 2 , z 3 , a 1 , \)L(z 3 , Zq, a 2 , A). (4.2) 

The left-hand side of the scalar equation (|3.3j ) also depends on the parameters: 

$(z ,zi,z 2 ,z 3 ,ai,a 2 ) = 0. (4.3) 

We shall require two additional features of the equation (|4.3| ). The first one is quite natural and 
is satisfied by all known examples: the function $ is symmetric with respect to the cyclic shift 
(zq, Zx, z 2 , z 3 ) i— ► (z\, z 2 , z 3 , zq) accompanied by interchanging the parameters a\ *-* a 2 . The second 
feature might seem not natural at all, at least on the first sight, but also holds, for some mysterious 
reasons, for the majority of (if not all) the known examples. Namely, suppose that the equation 
(|4.3| ) may be equivalently rewritten in the form 

f(zo, z 1 ,a 1 ) - f(zo, z 3 , a 2 ) = g(z , z 2 , ai,a 2 ). (4.4) 

We call the equation ( [4.4| ) the three leg form of the equation (|4.3| ). Indeed, its three terms may 
be associated with three "legs" (vo,v±), (vq,v 3 ), and (vq,v 2 ). Notice that while the first two "legs" 
belong to E(V), the third one does not; instead, it belongs to E{Q*) (recall that our quad-graph 
T> comes as a double constructed from G, G*)- Now consider n faces from F(T>) having common 
"black" vertex v G V(G) C V(V) (cf. Fig.|, where n = 5). Summing up the equations (|4.4| ) for 
these n faces, i.e. 

f(zo,Z2k-i,<*k) ~ f(zo,z 2 k+i,a k +i) = g(z ,z 2k ,ak,a k +i), (4.5) 
we end up with the equation which includes only the fields z in the "black" vertices: 

n 

^5(^o,22fc,afc,aA:+i) = 0. (4.6) 

k=l 

The latter formula may be interpreted as a system on the graph G- there is one equation per vertex 
vq G V(G), and this equation relates the fields on the star of vq consisting of the edges from E(G) 
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Figure 5: Faces from F(T>) around the 
vertex vq. 




Figure 6: The star (in Q) of the vertex 
Co- 



incident to vq , see Fig. ||. The parameters are then viewed as assigned to the corners of the faces 
of Q. Observe that Definition ^ of a labelling of E(T>) may be reformulated in terms of Q alone as 
follows: the corner parameters of two faces adjacent to the edge (t>o,f2fc) G E{Q) at the vertex vq 
are the same as at the vertex v 2k (the ordering of the parameters at vq when looking towards V2k 
is the same as the ordering of the parameters at V2k when looking towards vo, cf. Fig.^J). 

The collection of equations (^6) constitutes the discrete Toda type system on Q. Obviously, 
the nature of this system is formally somewhat different from the general construction of Sect. |2| 
(the former relates fields on stars, while the latter relates fields on closed contours). However, our 
construction suggests the following type of a zero curvature representation for ( [4.6|) . The wave 
function ^ : V(G*) i— * G[X] is defined on the "white" vertices V2k-i (in the notations of Fig. ||), i.e. 
on the vertices of the dual graph. The transition matrix along the edge e* = (v2k-i, ^2fc+i) £ E(Q*) 
of the dual graph is equal to 

L(z 2 k+i, Z2k,a>k, X)L(z 2 k-Z2k-i,ak+i, A) = L(z 2 k+i, z , otk+i, X)L(z , z 2 k-i, a>k, A). 

If one prefers to consider everything in terms of Q alone, one can think of as defined on the faces of 
Q, and of the above matrix as corresponding to the transition across the edge e = (vo, V2k) £ E(Q). 
One would like this matrix to depend only on the fields zq, Z2k- In all our examples this can and 
will be achieved by a certain gauge transformation ^f(vk) > A(vk)^{vk)- 

Next, notice that in the above reasoning the roles of Q and Q* may be interchanged. Hence, 
one can derive a Toda type system on Q* as well. The distribution of the corner parameters on Q 
and Q* is illustrated on Fig.[?]. 

We arrive at the following statement. 



Proposition 3 a) Let the system (4-j) on a quad-graph T> admit a three-leg form. Any solution 
z : V(V) i ^ C of this system, restricted to V(Q), resp. to V{Q*), delivers a solution of the discrete 
Toda type system (4-0 on Q, resp. on Q* , with the corner parameters given by the labelling of 
E{V). ~ ' 

b) Conversely, given a solution z : V{Q) i— > C of the discrete Toda type system on Q, there 
exists a one-parameter family of its extensions to a solution of the system ( f^. jj j on T>. 

Proof. The first statement is already proved. As for the second one, suppose we know the solution 
of the Toda type system on V(G), i-e. in all "black" vertices. Fixing the value of the field z in 
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Figure 7: Corner parameters of the dual graph 



one "white" vertex arbitrarily, we can extend it succesively to all "white" vertices in virtue of the 
equation ( [4.3D . The consistency of this procedure is assured by ( |4.4| ) . I 

Remark 4 Proposition || implies that from a solution of the Toda type system on Q one gets a 
one-parameter family of solutions of the Toda type system on Q* . This may be regarded as a sort 
of the Bdcklund transformation, especially when the graphs Q and Q* are isomorphic, which is 
the case, e.g., for the regular square lattice. 



5 Cross— ratio system and additive rational Toda system 
5.1 Definition 

The cross-ratio system is one of the simplest and at the same time the most fundamental and 
important integrable systems on quad-graphs. 

Definition 5 Given a labelling of E(D), the cross ratio system on T> is the collection of equa- 
tions for the function z : V(T>) i— > C, one equation per face of T> , which read, in notations as on 
Fig.$ and z k = z(v k ): 

, s def (z - z 1 )(z 2 - z 3 ) ai 

q(z ,z 1 ,z 2 ,z 3 ) = r = — . (5.1) 

{zi - z 2 ){z 3 - z ) a 2 



Several simple remarks are here in order. Although (|5.1| ) is a well-defined equation for each face 
of T>, the cross-ratio itself is not a function on F(T>), since 

q{z ,z 1 ,z 2 ,z 3 ) 



q(zi,z 2 ,z 3 .z ) 



On the other hand, F(T>) is in a one-to-one correspondence with E(Q), as well as with E(Q*). It 
is sometimes convenient to consider the cross-ratio as a function on E{Q) or on E(Q*), setting (see 
Fig.§ 

q(e) = q(z ,z 1 ,z 2 ,z 3 ), where e = (v , v 2 ) 6 E(G), (5.2) 
g(e*) = q(z 1 ,z 2 ,z 3 .z ), where e* = (v%, v 3 ) e E(Q*). (5.3) 

Sometimes a sort of inverse problem is of interest: given a function z : V{T>) i— * C, can it be 
considered as a solution of the cross-ratio system for some labelling? The answer follows from a 
simple lemma. 
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Lemma 6 Let, as usual, T> be a cellular decomposition of a topological disc (open or closed), and 
therefore a double for some pair of dual cellular decompositions Q , Q* . Let Q : E(Q) \—* C be some 
function, and extend it to E{Q*) by the formula Q(e*) = 1/Q(e). Then the necessary and sufficient 
condition for the existence of a labelling a : E(D) i— > C such that, in the notations of Fig. 

Q(e) = QK« 2 ) = — & 0(e*) = Q(«i,«3) = — , (5.4) 

Ct 2 0L\ 

is given (in the notations of Fig. ^) by the equations 

Yl Q{t) = 1 for all internal v G V(G), (5.5) 

eGstar(i>o) 

and 

Y\ Q(e*) = 1 for all internal vi G V(G*) (5.6) 

c* Sstar(t>i) 

Proof. The necessity is obvious. To prove sufficiency, we construct a by assigning an arbitrary 
value (say, a = 1) to some edge from E(T>), and then extending it successively using either of the 



equations (5.4) and the definition of labelling. The conditions ( |5.5| ) and ( p^q ) assure the consistency 
of this procedure. I 

Corollary 7 Let z : V(T>) i— > C be some map. The ncessary and sufficient condition for the 
existence of a labelling a : E(T>) C such that z is a solution of the corresponding cross-ratio 
system is given ( in the notations of Fig. [^) by the equation 

n 

Y\ l( Z 0, Z2k~l,Z 2 k, Z 2k+ l) = 1 
fe=l 

for all internal vertices vq G V(T>). 

5.2 Zero curvature representation 

Proposition 8 The cross-ratio system on T> admits a zero curvature representation in the sense 
of Sect. [J ; with the transition matrices for e = (vq ,v\): 

( 1 z - Z\ \ 

L(e, A) = L( Zl ,z , a, A) = (1 - Aa)" 1/2 



\ z - Z\ 



(5.7) 



Here a = a(e), and, as usual, z/% = z(vk). 



Proof. Verifying Q4.1| ) and ( f4.2| ) is a matter of a simple computation. Notice that the entry 12 
of the matrix equation (fl]^) gives an identity, the entries 11 and 22 give the equation ( |5.1| ) as it 
stands, while the entry 21 gives: 

«1 «2 «2 Oil 

H — 1 , (o-oj 



z - Z\ z\- z 2 z - z 3 z 3 - z 2 



which is a consequence (5.1). I 

The zero-curvature representation of Proposition ||| was introduced in [NC|, of course only in 



the case of the square lattice and constant labelling. Cf. also a gauge equivalent form in |BP1 
where a more general situation of discrete isothermic surfaces in M 3 is considered (in this context 
the fields z take values in H - the ring of quaternions). 
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5.3 Zero curvature representation from three-dimensional consistency 

Now we demostrate how to derive the zero curvature representation of Proposition ||| starting from 
an additional deep property of the cross-ratio system. To this end we extend the planar quad- 
graph T> into the third dimension. Formally speaking, we consider the second copy T>' of T> and 
add edges connecting each vertex v G V(T>) with its copy v' £ V(V). On this way we obtain a 
"three-dimensional quad-graph" D, whose set of vertices is 



whose set of edges is 



V(D) = V(V)UV(V), 
E(D) = E(V) U E(V) U {(«, »'):»£ V(V)}, 



and whose set of faces is 

F(D) = F(V) U F(V) U {(^o^i^i^o) : v o,vi G V(V)}. 
Elementary building blocks of D are "cubes" (vq, v\, V2, v%, v' , v[, v' 2 , v' 3 ), as shown on Fig.|8|. Clearly, 



a 



(i 
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Figure 8: Elementary cube of V U T>' 



we still can consistently subdivide the vertices of D into "black" and "white" one, assigning to each 
v' G V(V) the colour opposite to the colour of its counterpart v G V{T>). 

Next, we define a labelling on E(D) in the following way: each edge (v' ,v[) € E(D') carries 
the same label as its counterpart (^0,^1) G E(V), while all "vertical" edges (v,v ! ) carry one and 
the same label \i. Now, the fundamental property of the cross-ratio system mentioned above is the 
three dimensional consistency. This should be understood as follows. Consider an elementary 
cube of D, as on Fig.S. Suppose that the values of the field z are given at the vertex vq and at its 



three neighbors v\, v%, and v' . Then the cross-ratio equation (5.1) uniquely determines the values 
of z at V2, v[, and v%. After that the cross-ratio equation delivers three d priori different values 
for the value of the field z at the vertex v' 2 , coming from the faces (vi,V2,v 2 ,v[), (v2,V3,v' 3 ,v' 2 ), 
and (v' , v[, v' 2 , v' 3 ), respectively. The three-dimensional consistency means exactly that these three 
values for z{v' 2 ) actually coincide. The verification of this claim consists of a straightforward 
computation. Because of the importance of this property we formulate it in a separate proposition. 



Proposition 9 The cross-ratio system is consistent on the three-dimensional quad-graph D. 

Now we can derive a zero curvature representation for the cross-ratio equation on T>, using the 
arbitraryness of the label [i assigned to the "vertical" edges (v,v') of D. For this aim, consider 
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the cross-ratio equation on the "vertical" face (vq, v\, v[, v' ). Denote A = \x , then the equation 
reads: 



(z[ - z' q )(zq - Zl) 
(z' -z )(zi -z[) 



Xa. 



(5.9) 



This gives z[ as a fractional-linear (Mobius) transformation^] of z' with the coefficients depending 
on zq, zi, and the (arbitrary) parameter A: 



where 



z'i = L(zi,zo,a,X)[z' ], 



/ ^ Xazi XozqZi \ 



L(z 1 ,z ,a,X) 



z - Zl 
Xa 



z - Zl 



Z - Zi 

Xazo 

Zq - Zl 



1 + 



Xa 



zq - Zi \ 1 



-Zq 



(1 - Xa)I + 



Xa ( zo 



-zi 



(5.10) 



(5.11) 



In principle, the matrices L(zi, Zq, a, A), after multiplying by (1 — Aa)" 1 / 2 (which is nothing but 
normalizing the determinants) are the sought after transition matrices of the discrete flat connection 
on T>, as assured by Proposition ^. To recover the transition matrices from Proposition || we 
perform a gauge transformation defined on the level of the wave functions : V(T>) \— > G as 



i-> A(z(v))W(v), 



where 



A(z) 



1 z 
1 



(5.12) 



(5.13) 



On the level of transition matrices this results in 



L(zi,z ,a, A) i-> L(zi,z ,a, A) = yl z , «, A)/l(zo) 



/ 1 

Xa 



zq - zi 



\ Zq- Zl 



1 



/ 



Notice that the gauge transformation ( 5.12[) , ( |5.13 ) may be interpreted in terms of the variables 
z(v') as a shift of each z(v') by z(v), so that the interpretation of the matrix L(zi, Zq, a, A) is as 
follows: 

z[- zi = L(z 1 ,z , a, X)[z' - zo]. 



We use the following matrix notation for the action of Mobius transformations: 



az + b 
cz + d 



L[z], where L 



a b 
c d 
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From our point of view, the property of the three-dimensional consistency lies at the heart 
of the two-dimensional integrability. Not only does it allow us to derive the spectral-parameter 
dependent zero-curvature representation for the two-dimensional discrete system. We can consider 
the function z : V(T>') h C as a Bdcklund transformation of the solution z : V(T>) C. 
Moreover, if T> is generated by a square lattice, then, as well-known, by refining this lattice in 
one or two directions, we can obtain certain integrable differential-difference or partial differential 
equations, the latter being the famous Schwarzian Korteweg-de Vries equation 

3 z xx 

Zt — Z xxx i 
Z Z x 



sec 



NC|. So, the three-dimensional consistency delivers all these equations and their Backlund 
transformations in one construction. For instance, if in ( |5.9| ) zq is a solution of the Schwarzian 
KdV, z\ its Backlund transformation, and z' Q , z[ correspond to the e-shift of these solutions in the 



^-direction, where A = e , then in the limit e — > one recovers the well-known formula | Al | 

(zo)x(zi) x = ol~ 1 {zq - zx) 2 . 
Similarly, when considered on a single face of the diamond, as on Fig.^j, the cross-ratio equation 



(5.1) expresses nothing but the commutativity of the Bianchi's diagram (Fig.|4|) for the Backlund 



transformations of the Schwarzian KdV. Indeed, multiplying the equations 

(zo)x(zi)x = a^ 1 (z - zi) 2 , {z 3 ) x (z 2 ) x = a{ 1 {z 3 - z 2 ) 2 , 
and dividing by two similar equations corresponding to the edges (vq,v 3 ) and (^1,^2) (with the 



parameter a 2 ), we arrive at the (squared) equation (5.1). So, one can say that the equation ( 



on the cube Fig. |8| contains everything about the Schwarzian KdV and its Backlund transformations. 
5.4 Additive rational Toda system 

The further remarkable feature of the cross-ratio equation is the existence of the three-leg form. 
To derive it, notice that (|5.1| ) is equivalent to 



or to 



z\ - z 2 z 3 - z 2 z - z 2 z - z 2 

ct\ a 2 = U 44> a\ a 2 = ct\ — a 2 , 

zo ~ z\ z - z 3 z -z\ zq- z 3 



qi a 2 ai - a 2 

(5.14) 



Z - Z\ Zq - Z 3 Zq - Z 2 

which obviously is of the form (f4,4|). 

According to the construction of Sect.Q, the corresponding system of the Toda type reads (in 
the notations of Fig. ||) : 

n 

y a k - a k+1 = Q 

It will be called the additive rational Toda system on Q. For the cross-ratio system on T> and 
additive rational Toda systems on Q, Q* there holds Proposition |3|. 

Remark 10 Notice that the Toda system on Q remains unchanged upon the simultaneous additive 
shift of all corner parameters a k 1— ► Ofc + a. However, such a shift changes the cross-ratio system 
essentially. Hence, Proposition || yields also a nontrivial transformation between the solutions of 
two different cross-ratio systems on T> whose labellings differ by a common additive shift. 
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Next, we derive a zero curvature representation for the additive rational Toda system from such 
a representation for the cross-ratio system. Consider the additive rational Toda system on the 
"black" vertices from V(Q), i.e. on v 2 k m the notations of Fig. ||. The transition matrix across the 
edge e fc = {v ,v 2k ) G E(Q), i.e. along the edge z* k = {v 2k -i,v 2k+ i) G E(Q*), is given by 

£(z 2k , Zo, A) = L(z 2k+1 ,z , a k+ i, X)L(zo, z 2k _ x ,a k , A). (5.16) 

We use here the matrices fl5.1l|) , gauge equivalent to (|5.7| ), because of the remarkable fact underlined 
by the notation for the matrices on the left-hand side of the last equation: they actually depend 
only on zo,z 2k , but not on z 2 k-i, z 2k+ \. (This propery would fail, if we would use the matrices 

(0)0 

Proposition 11 



ak+iz 2 k - akZQ [ak - ak+i)z z 2 k 
C(z 2k , zq, A) = I + —— | j (5.17) 

«fc+i - OL k OL k z 2k - a-k+izo 

= (l-Ao fc )/ + A»f Z f)( 1 V (5.18) 
zo ~ Z 2 k V 1 / V ~ z o J 

= (1-Aa fe+1 )/ + A^±i^f * )( 1 Y. (5.19) 

ZO - Z 2k V 1 / V ~ Z 2fc / 

Proof. Use the second formula in ( 5.1 lj ) for L(zq, z 2 k—ii a k-> A) and the first formula in ( |5.11| ) for 
L(z 2 k+i, Zo, ojfe+i, A). This leads, after a remarkable cancellation of terms quadratic in A, to the 
following expression for the matrix ( [5.16| ): 



£(z 2k , z , A) = (1 - \a k )I + A£ 

where 



1 

-zo 



^ _ afc+i ( z 2k+ i \ a>k ( z 2k -i 



zo - z 2k +i \ 1 / z - z 2 k-i V 1 
A simple calculation based on the formula ( |5.14 ), written now as 

Qfc+l «fc _ Qfc+l - «fc 

— ^2A:+l ^0 — Z 2 k-1 Zo ~ Z 2 k 

shows that 

«fc+i -otk I z 2k 



zo - z 2k V 1 
This proves d5.18|) . I 

The matrix ( |5.17 ) coincides with the transition matrix for the additive rational Toda lattice 
found by V. Adler in [ ASf] . 



13 



5.5 Duality 

Next, we discuss a duality transformation for the cross-ratio system. This comes into play, if we 
rewrite the cross-ratio equation (|5.8| ) in the form 



_^ + ^^ + ^i_ + ^^ = . (5.20) 

Zq - Zl Z\- Z 2 Z 2 -2 3 Z 3 - Zq 

This guarantees (locally, and therefore in the case when D is a cellular decomposition of a topological 
disc) that one can introduce the function Z : V(T>) i— » C such that 

Zi — Zq = — , Z 2 — Z\ = — , Z3 — Z 2 = — , Zq — Z3 = — . (5.21) 

ZQ ~ Zl Zl - Z 2 Z 2 ~ Z3 Z3 - Zq 

In other words, for an arbitrary e = (t>i,t>2) € E(T>), one sets 

Z(v 2 )-Z(vi)= " (e) (5.22) 

z(Vl) - Z(V2) 



Proposition 12 For any solution z : V(T>) i—* C of the cross-ratio system on T> , the formula 



5.2% ) correctly defines a new (dual) solution of the same system. 



Proof. The above argument shows that the formula ( |5.22 ) correctly defines some function on 



V(T>). It can be easily seen that this function also solves the cross-ratio system. I 

Combining Propositions ||| and 12, we see that to any solution of the additive rational Toda 



system on Q there corresponds a one-parameter family of solutions of the additive rational Toda 
system on Q* , and any such pair of corresponding solutions yields (almost uniquely - up to an 
additive constant) dual solutions for the additive rational Toda systems on Q and Q*. But actually 
this duality relation between Toda systems is much more stiff. Indeed, adding the first and the 
fourth equations in ( |5.21 ) and using ( 5.14| ), we arrive at the formula 



Zi-Z 3 = (5.23) 

Z2 ~ Zq 

Therefore, a solution of the additive rational Toda system on Q determines the dual solution of the 
additive rational Toda system on Q* almost uniquely (up to an additive constant). 

It is also instructive to regard the function w : V{T>) ^ C which coincides with z on V(Q) 
and with Z on V(Q*). Then formula ( |5,23j ) may be considered as a system of equations on the 
quad-graph, reading on the face (vo,vi,V2,vs) as 

a 2~a>i . . 

wi — u>3 = . (5.24) 

W2 - Wq 



Interestingly (but not very surprizingly) , the equations ( [5.24 ) also constitute an integrable system 



on the quad-graph. In the case when T> is generated by a square lattice and the labelling is constant, 
this system is known under the name of the "discrete KdV equation" [ jHl| |, [ NC| . Actually, it 
appeared even earlier in the numerical analysis as the "e-algorithm" |W| . The transition matrices 
for the equation ( 5.24| ) are given by 



wq —wqWi + a — // 

L(e,/i) = L(w\,wq,ol,h) = (p - a)~ 1/2 ( |. (5.25) 

1 — Wi 
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This transition matrix may be derived in exactly the same way as for the cross-ratio system, based 
on the three-dimensional consistency. Namely, this is the matrix of the Mobius transformation 



-» v[, see Fig. ^: 

(w\ - w' )(w[ - w ) + a — /j, = 



wqWq — wqWi + a — fj, 



L(wi, w , a, fj,)[w' ]. 



■j -wi 

Again, the three-dimensional consistency is the central principle for the whole integrability theory 
of the discrete, differential-difference and partial differential equations related to the (potential) 
Korteweg-de Vries equation 

„2 



W t 



w 3 



as well as their Backlund transformations. In particular, the equation fl5.24|) expresses nothing 
but the commutativity of the Bianchi's diagram (Fig.^J) for the Backlund transformations of the 
potential KdV. Indeed, recall that the Backlund transformations corresponding to the edges (vq, V\) 
and (^3,^2) of the diagram Fig.|3] are defined by the equations 

{wq) x + {Wi) x = Oi - W ) 2 + Oil, {W 3 ) x + (W 2 ) x = (W 2 ~ W 3 f + Qi. 

Adding these two equations and subtracting two similar equations corresponding to the edges 
(^0,^3) and (vi,V2) (with the parameter 02), we arrive at 

(u> 2 - w )(wi - W3) + ai - a 2 = 0, 



which is nothing but ( 5.24 ). 



Notice that equation ( 5.24 ) is already in the three-leg form, and restrictions of its solutions to 
V(G) and to V(Q*) are solutions of the rational additive Toda system on the corresponding graphs. 
It is easy to check that the pairwise products of the matrices ( |5.25| ), like in ( 5.16| ), lead to the same 
transition matrices (|5.17 ) for the additive rational Toda system as before. 

The relations between various systems discussed in the present section are summarized on Fig. || 



Toda on Q 



discr. KdV on T> 



cross-ratio on T> 





duality 


duality 





Toda on Q* 



discr. KdV on T> 



dual Toda on Q* 



dual cross-ratio on V 



dual Toda on Q 



Figure 9: Relations between different systems 



The thick lines here are for the restriction maps, the thin ones - for the duality maps. The 
duality for the cross-ratio system is described by the formula (5.22). Its transcription in terms of 
the variables w and the dual ones W reads: 

a(e) 



W(V2) — w(v\) 



W{y\) — w{v2) 



(5.26) 
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It is not difficult to prove the statement similar to Proposition |T^: For any solution w : V(T>) \— > C 
of the discrete KdV system on D, the formula ( 5.26 ) correctly defines a one-parameter family of 
(dual) solutions W : V(T>) i— > C of the same system. 



5.6 Cross— ratio system and circle patterns 

We close the discussion of the cross-ratio system on an arbitrary quad-graph T> by giving its 
geometric interpretation. For this aim, we consider circle patterns with arbitrary combinatorics. 

Definition 13 Let Q be an arbitrary cellular decomposition. We say that a map z : V(Q) i— > C 
defines a circle pattern with the combinatorics of Q, if the following condition is satisfied. 
Let f G F(Q) be an arbitrary face of Q, and let vi,v%, . . . ,v n be its consecutive vertices. Then the 
points z(v\), z(v2), • • • , z(v n ) G C lie on a circle, and their circular order is just the listed one. We 
denote this circle by C(f), thus putting it into a correspondence with the face f, or, equivalently, 
with the respective vertex of the dual decomposition Q* . 

As a consequence of this condition, if two faces /i,/2 G F(Q) have a common edge (vx^v-f), then 
the circles C(/i) and Ci^f-f) intersect in the points z(vi), z{v2)- In other words, the edges from 
E{Q) correspond to pairs of neighboring (intersecting) circles of the pattern. Similarly, if sev- 
eral faces /i,/2,... , fm ^ F(Q) meet in one point vq G V(G), then the corresponding circles 
C(/i), C(ff), • • • , C(f m ) also have a common intersection point z(vq). 



Given a circle pattern with the combinatorics of Q, we can extend the function z to the vertices 
of the dual graph, simply setting 




Figure 10: Circle pattern 



z(f) = Euclidean center of the circle C(f), f G F{Q) ~ V(G*). 



A slightly more general extension is achieved by setting 



z{f) = conformal center of the circle C(/), / G F(Q) ~ V(Q*). 
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Recall that the conformal center of a circle C is a reflection of some fixed point G C in C. In 
particular, if P m = oo, then the conformal center coincides with the Euclidean one. 

Thus, after either extension, the map z is defined on all of V{T>) = V(Q) U V(Q*), where V 
is the double of Q. Consider a face of the double. Clearly, its vertices vq,v\,V2,v^ correspond to 
the intersection points and the (conformal) centers of two neighboring circles Co, C± of the pattern. 



Let zq, Z2 be the intersection points, and let z\, z% be the centers of the circles. See Fig. 11 




Figure 11: Two intersecting circles 



Lemma 14 If the intersection angle of Cq,C\ is equal to cj), then 

q(z , Zi,z 2 , z 3 ) = exp(2icj)). (5.27) 

Proof. The claim is obvious for the Euclidean centers. These can be mapped onto conformal 
centers by an appropriate Mobius transformation which preserves the cross-ratios. I 

Proposition 15 Let Q be a cellular decomposition of a topological disc (open or closed), and con- 
sider a circle pattern with the combinatorics of Q. Let z : V(T>) i— > C come from a circle pattern, 
i.e. {z(v) : v G V{G)} consists of intersection points of the circles, and {z(v) : v G V{Q*)} consists 
of their conformal centers. Suppose that the intersection angles of the circles satisfy the following 
condition: 

• For each circle of the pattern the sum of its intersection angles 4>k, k = 1,2, ... ,m with all 
neighboring circles of the pattern satisfies 2 X^fcLi 0fc = (mod 2tt) . 

Then there exists a labelling a : E(D) i— > S 1 = {6 G C : \9\ = 1} depending only on the intersection 
angles, such that z is a solution of the corresponding cross-ratio system. 

Proof. According to Corollary M, the necessary and sufficient condition for the claim of the 
proposition to hold is given by (|5.5| ), (|5.6|). The second of these two formulas is exactly the 
condition formulated in the proposition. The first one is equivalent to the following condition: for 
an arbitrary common intersection point of n circles the sum of their pairwise intersection angles 
01,... , 4> n satisfies the similar condition 2j^" =1 0j = (mod 2tt). This condition is satisfied 
automatically for geometrical reasons (its left-hand side is equal to 2ir). I 

Corollary 16 Under the condition of the previous proposition, the restriction of the function z to 
V{Q) (i.e. the intersection points of the circles) satisfies the additive rational Toda system on Q, 
while the restriction of the function z to V{Q*) (i.e. the conformal centers of the circles) satisfies 
the type (A) Toda system on Q* . 
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6 Shifted ratio system and multiplicative rational Toda system 



We now define a system on a quad-graph T> which can be considered as a deformation of the 
cross-ratio system. 

Definition 17 Given a labelling of E(T>), the shifted ratio system on T> is the collection of 
equations for the function z : V(T>) i— ► C, one equation per face of T>, which read, in notations as 
on Fig. and = z(vh): 



(zq - zi + ai){z 2 - z 3 + a{) 



a i 

0' 2 



(6.1) 



Ol 
«2 



(6.2) 



(zi - z 2 - a 2 ){z 3 - z - a 2 ) 

or, equivalently (as one checks directly), 

(z - z 1 - ai)(z 2 - z 3 - ax) 
(zi - z 2 + a 2 )(z 3 - zq + a 2 ) 

The cross-ratio system of the previos section may be obtained from the shifted ratio system by the 
following limit: replace in (6J.) by eat and then send e — > 0. The theory of the shifted ratio 
system is parallel to that of the cross-ratio system, therefore our presentation here will be more 
tense. 

Proposition 18 The shifted ratio system on T> admits a zero curvature representation in the sense 
of Sect. p| with the transition matrices for e = (vq,vi): 



( 



L(e, A) = L(zi, zq, a, A) = (1 — Aa) 



-1/2 



1 

Aa 



zq 



Z\ + a \ 
z\ + a 



(6.3) 



\ Zq — Z\ — a Zq — Zi — a / 



Here a = a(e), and z^ = z{vk). 

Proof. As in the proof of Proposition ||, we have to verify the equations ( |4.1| ) and (|4.2| ). This is 
done again by a direct computation. In particular, one checks that the entry 11 and the A-term of 
the entry 22 of the matrix equation 

L(z 2 , z 1} a 2 , X)L(zx,zo, «i, A) = L(z 2 , z 3 , at, X)L(z 3 ,z , a 2 , A). 

are equivalent to the equations (|6.2| ) and (|6.1| ), respectively. The A°-term of the entry 22 reads: 
Z\ — z 2 + a 2 z - zi+ ai z 3 - z 2 + u\ z - z 3 + a 2 



Zi — z 2 — a 2 zq - z\- a\ z 3 - z 2 - a\ z - z 3 - a 2 
which is a plain consequence of fl6.1|) , (|6.2[) . Further, the entry 12 reads: 



(6.4) 



z - z\ + a\ + (zi - z 2 + a 2 ) 



Zq — Zi + ai 

Zq- Zi- ai 



zq - z 3 + a 2 + (z 3 - z 2 + ai) 



Zq - Z 3 + a 2 
Zq- Z 3 -a 2 



which is equivalent to 



zq — z\ + a\ zq — z 3 - a 2 _ zq — z 2 + a\- a 2 
Zq - Zi - a\ Zq - z 3 + a 2 zq - z 2 - a\+ a 2 



(6.5) 
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The latter equation (after taking logarithm) is nothing but the three leg form of the shifted ratio 
equation. Indeed, straightforward algebraic manipulations show that it is equivalent to ( |6.1| ) (or to 
|)). Finally, consider the entry 21 of the above matrix equation: 

CK2 OL\ Z\ - Z 2 + «2 _ "1 «2 ^3 ~ z 2 + «1 



Zl — ^2 — "2 Zq — Zi — Oi Zi - Z 2 - «2 ^3 - ^2 - "l ^0 ~ z 3 ~ «2 ^3 ~ z 2 ~ Oil 

With the help of ( |6.2| ) this may be rewritten as 

a 2 ( l z 2 - z 3 - qi \ _ ai / zi - z 2 + a 2 



Zl — ^2 — «2 V -^3 - -20 + "2 / -^3 - -22 - "1 V 2 — &l ~ a l 



and this, in turn, is rewritten with the help of ( |6.1[) as (6.5). I 

The shifted ratio system, as well as its generalization from the next Section, appeared in | NQC 
(of course, only in the simplest possible situation: regular square lattice, labelling constant along 
each of the two lattice directions). They denoted this equation as "the discrete Krichever-Novikov 



system". The zero curvature representation was given later in [NC|. Actually, this zero curvature 
representation can be derived along the same lines as for the cross-ratio system, based on the 
three dimensional consistency, which remains to hold true also in the present situation. 

Proposition 19 The shifted ratio system is consistent on the three-dimensional quad-graph D. 

Proof is again a matter of straightforward computations. I 

To derive a zero curvature representation from this fact, we proceed as in the previous section: 
assign the label [i to all "vertical" edges (v,v f ) of D, and consider the shifted ratio equation on the 
"vertical" face (vq ,Vi,v[,v' q ): 

(z[ - z' + a)(z — zx + a) _ a 
{z' - z - y)(z\ - z[ - n) fi' 

Denote A = fj, , and solve for z' x in terms of z' : 

z[ = L(z 1 ,z ,a,X)[z' Q ], 

where 

' zq — z\ + Xa X~ 1 a — a 2 - 



L(z 1 ,z ,a,X) 



Xa Zq — Z\ 




(Recall that the limit to the formulas of the previous section is the following one: a ^ ea, /i h- > e/x, 
so that A i — > e _1 A, and then e — > 0.) To recover the matrices of Proposition perform the gauge 
transformation 



h-» C(z(v),X)tf{v), (6.6) 



where 



1 z-X- 1 \ 






r ) 


1 J 


H 


V o 





C(z,X)= = \A(z). (6.7) 
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An easy computation shows: 



zq — z\ — a (zq — zi) 2 — a 2 

L{z 1 ,z ,a,X) = C~ 1 {zi,X)L{zi,z ,a, X)C(z , A) 

Xa zq — z\ + a 

which coincides with (|6.3| ) up to a scalar factor (inessential when considering action of L as a 
Mobius transformation) . 

Since the shifted ratio equation possesses the three-leg form, there holds the analog of Propo- 
sition ||. The corresponding systems of the Toda type on Q, Q* will be called the multiplicative 
rational Toda systems. Their equations (in notations of Figs. read: 

n i 

TT ZQ ~ Z 2 k -Oik + =1 

C i z o - z 2 k + at k - a k +i 

k=l 

As in the previous section, we can find the transition matrix across the edge = (yo,V2k) £ E(Q), 
i.e. along the edge z* k = {v2k-i-i v 2k+\) £ E{Q*). This time we start with the "untilded" matrices 

£fc(A) = L(z 2 k+i, zo, afc+i, A)L(z , ^fc-i, «fc, A). (6.9) 



Proposition 20 Under the inverse gauge transformation \5.1i ), (5.11), and up to the scalar de- 
terminant normalizing factor, the matrix turns into 

( , x Zq + Z 2 k ~ Ot k - Oik+1 \ 
I 1 {ak+i ~ OL k ) ; \ 



£{z 2 k,z ,X) = A(z2k+i)Ck{X)A 1 (z 2k -i) 







zq 



zo - Z2k + otk - ak+i 

Z2k - Oik + Q-k+1 



+ 



+ - 



A 



\ zq - z 2 k + a k - a k+ i 

ot-k+\Z2k - oikZQ {ak - ak+i){z z 2 k - afcQ!fe+i) 



ZQ - Z2k + OLk - Otk+1 



(6.10) 



ttfc+l - OLk 



OtkZ2k - OLk+\ZQ 



Proof. The calculation of the A°-component of the matrix ( |6.10 ) is immediate, one has only to 
use the three-leg equation (6.5), i.e. 



ZQ - Z2k+\ + ak+1 _ Zq - Z 2 k~l ~ Qk _ Zq - Z 2k ~ «fc + a k +l 

zq - z 2 k+i - ak+i zq - z 2 k-i + a k Zq - z 2 k +ak- «fc+i 



(6.11) 



The calculation of the A-component is also straightforward but somewhat more involved. The main 
ingredient of this calculation is obtaining the following expression for the entry 21 (which is used 
actually in all four entries): 



zo — Z2k+i + ak+i 



+ 



ak+i 



ak+i - ak 



Z2k-i - zo- ak z - z 2 k+i - ak+i zo - z 2 k+i - ak+i z 
This is done as follows: one starts with the shifted ratio equation 



Z2k + a k - a k+ i 



(6.12) 



a k 



z 2k — ggfc— 1 — ak+1 
Z2k-1 ~ Zq + a k 



ak+i 



z 2k+l — z 2k + ak 

zo — Z2k+i — ak+i 



0. 
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This is rewritten as 



zq - z 2k - a k + a k+ i z 
OLk ; ; — — V oik+i 



z 2 k + otk - at k+ i 



Oik+l - Oik- 



z 2 k-\ - zo + a k Zq — z 2 k+i - afc+i 

But the latter equation is equivalent to ( |6.12| ) due to the three-leg equation ( |6.11 ), 



The transition matrix for the type (C) Toda system found in [A3| coincides with (6.10) up to 
a similarity transformation with a constant matrix I ^ 

Finally, we briefly discuss the dual to the shifted ratio system. It appears, if one notices that 
the equation (6.4) yields the existence of the function Z : V(T>) t— ► C such that on each face of T> 
there holds: 



ry ry 1 , Z Q ~ Z 1 + Oil 

I Zq — Z\ — OL\ 



7 7 1 , z x - z 2 + a 2 
Z 2 - Zi = - log . 

2 z\ - z 2 - a 2 



ry ry 1 , Z 2 - Z 3 + Qj 

Z 3 -Z 2 = - log . 

2 z 2 - z 3 - ai 



ry ry 1 , ^3 ~ Z + Q 2 

Z - Z 3 = - log . 

2 23 - Z - Oi 2 



Proposition 21 For any solution z : V(2?) i— »■ C o/ i/ie shifted ratio system on T>, the formula 

Z(v 2 )-Z{vi) = -\og—— — — (6.13) 

2 z(y 2 ) - z(yi) - a(t) 

(for any e = (v\,v 2 ) € E(T>)) correctly defines a function Z : V(T>) i— > C which is a solution of the 
hyperbolic cross ratio system 

sinh(Z - Zi) sinh(Z 2 - Z 3 ) _ oti 
sinh(Zi — Z 2 ) sinh(Z3 — Zq) ct 2 

Remark 22 The hyperbolic cross-ratio system is equivalent to the usual cross-ratio system 

(X - X 1 )(X 2 - X 3 ) _ at 
{X l - X 2 )(X 3 - X Q ) a 2 

upon the change of variables X& = exp(2Zfc). Restrictions of solutions of the hyperbolic cross-ratio 
system to V(Q), V{Q*) satisfy the additive hyperbolic Toda systems, whose equations are 

n 

^2 («k - «k+i)coth(2;o - z 2k ) = 0. (6.15) 
k=l 

This system is also equivalent to the additive rational Toda system 

n 

sr^ a k - Ot-k+l _ q 
x - x 2k 

under the change of variables x& = exp(2zfc). 
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7 Hyperbolic shifted ratio system and multiplicative hyperbolic 
To da system 

Finally, we turn to the analog of the cross-ratio system which is the most general one. 

Definition 23 Given a labelling of E(T>), the hyperbolic shifted ratio system on V is the 

collection of equations for the function z : V(T>) i— > C, one equation per face of T>, which read, in 
notations as on Fig. [| and Zk = z(vk): 

sinh(zo — Z\ + a\) sinh(;22 — zj, + ol\) sinh(2ai) (7 1) 



or, equivalently, 



sinh(,zi — Z2 — 0.2) sinh(2;3 — zq — 02) sinh(2a2) 
sinh(z — Z\ — at\) sinh(z2 — £3 — a{) sinh^ai) 



sinh(,zi — Z2 + 02) sinh(z3 — zq + 02) sinh(2a2) 



(7.2) 



The hyperbolic shifted ratio system may be considered as a deformation of the rational one, 
the corresponding limit from the former to the latter is achieved by setting replace z^ ^ ez^ and 
afc — > eak and then sending e — > 0. 

We shall systematically use in this section the notations Xk = exp(2zk), Ck = exp(2a/c). In these 
notations the equations ( |7.1[ ), ( |7.2|) take the following rational form: 

(CjXQ - xi)(cix 2 - x 3 ) _ 1 - cf 
(x\ - c 2 x 2 )(x 3 - c 2 x ) l-cl ' 

and 

(x - c 1 x 1 )(x 2 - cix 3 ) 1 — 



(C 2 X1 - X 2 )(C2X 3 - Xq) 1-C 



2 



(7.4) 



Proposition 24 The hyperbolic shifted ratio system on T> admits a zero curvature representation 
in the sense of Sect.^, with the transition matrices for e = (vq,vx): 



I 1 Xq-CXx \ 

i(,,A) = L(x 1 ,x , c ,A) = (l-A(l-c^ xa _^ 



(7.5) 



\ CXq — Xi CXq — X\ J 



Here c = c(z), and Xk = x(vk)- 

Proof is absolutely parallel to the proof of Proposition I 

Again, this statement may be derived from the three dimensional consistency of the shifted 
hyperbolic ratio system. 

Proposition 25 The hyperbolic shifted ratio system is consistent on the three-dimensional quad- 
graph D. 
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Proof - a direct computation. I 

The derivation of a zero curvature representation from this fact is not different from the previous 
two cases: we assign the label fi to all "vertical" edges (v,v') of D, set v = exp(2/i), and consider 
the hyperbolic shifted ratio equation on the "vertical" face (vq,vi,v'i,v'q): 



(cxq - xi)(cx[ - x' ) 



1 



(x\ — ux f - i )(x' — UXq) 1 — V 2 
Solving for x\ in terms of x' , we find: 

x[ = L(xi,x ,c,v)[x' ], 

where 



/ 1 

' CXq —X\ + 



L(xi,x ,c,v) 



Perform a gauge transformation 



where 



1/(1 



1-u 

„2\ 



-XqXi 



v- 



Xq — CX\ 



7^0 



i-> D{x(v),v)^(v), 



V 1 V l X 



D(x,v) 





o\ 


)-l 


V o i j 



(7.6) 



(7.7) 



L(xi, x , c,X) = D (xi,f)L(xi,x , c, j/)-D(xo, v) 



1 

It leads to the following form of the transition matrices: 

cxo — x\ (cxo — x\)(xo — CX\) 

A(l — (?) xo — cx\ 

where A = (1 — v 2 ) -1 . This coincides with ( [7.5[ ) up to an inessential scalar factor. 

Next, we turn to the derivation of the Toda type system related to the hyperbolic shifted ratio 
system. As usual, the three leg form of the latter system is fundamentally important in this 
respect. It is not difficult to check that ( |7. 1|) is equivalent to the following equation: 

sinh(zo — %\ + «i) sinh(z - z 3 - a 2 ) _ sinh(z - z 2 + ai - a 2 ) g , 
sinh(zo — z\ — a.\) sinh(zo — z% + a 2 ) sinh(zo — z 2 — oi\ + a 2 ) 

The rational form of the three-leg equation reads: 

{ciXq - xi) _ (x Q - c 2 x 3 ) _ {cixp - c 2 x 2 ) ^ ^ 

(xQ-ciXi) (c 2 x -x 3 ) (c 2 x -cix 2 )' 

The three-leg form of the equations immediately yields the analog of Proposition ||. The corre- 
sponding Toda type systems on Q, resp. Q* , read (in notations of Figs. 



n 



sinh(z - z 2k + a k - ctk+i) 
L sinh(z - z 2k - a k + a k+1 ) 



(7.10) 
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or in the rational form 

n 

l r CfcXp - c k+ \x 2k _ ^ ^ ^ 

k=i Ck+lX ° ~ Ck0 ° 2k 

They will be called multiplicative hyperbolic Toda systems. 

Turning to the zero curvature formulation of the multiplicative hyperbolic Toda system on G, 
we find, as in the previous section, the transition matrix across the edge c k = (vo,v 2 k) £ E(G), i.e. 
along the edge e* k = (v 2k -i,v 2 k+i) G E(Q*)\ 

Ck(X) = L(x 2k+ i,x , c fe+ i, X)L(x Q , x 2 k-i,ck, A). (7.12) 



Proposition 26 Under the inverse gauge transformation ^5.1 r k ), (5.1c), and up to the scalar de- 
terminant normalizing factor, the matrix ft7.13J turns into 



C(x 2k ,x ,X) = A(x 2k+1 )C k (X)A (x 2 k-i) 



^ (1 - A) + Xc k c k+ i- k ' " 



C k +\X 2k 



c k+1 x - c k x 2k 



(1-A)- 



-k+l 



" k )X[)X 2k 



Cfc+iXQ - C k X 2k 



-k+l 



C k+1 X - C k X 2k 



(1-A) 



c k x - C k+ lX 2k 
C k+ iX - c k x 2k 



+ Xc k c k+ i 



(7.13) 



Proof - by a direct calculation. The key formula in this calculation is: 



xq — C k+ \X 2k +l 



1 



(7.14) 



XO - C k X 2k -\ C k+ \XQ - x 2k+ i c k+1 x - X 2k+1 C k+ lX - c k x 2k 

which is proved similarly to fl6.12|) , based on the three-leg formula (|7.9| ). I 

Again, the transition matrices for the multiplicative hyperbolic Toda systems essentially equiv- 
alent to ( [T.13D were found in |A3| . 



The multiplicative hyperbolic Toda lattice can be given a geometrical interpretation. For this 



aim consider, as at the end of Sect. 5.6, a circle pattern with the combinatorics of a cellular 



decomposition G, and extend it by the Euclidean centers of the circles. So, the Euclidean radii of 
the circles r : V(G*) | — ► R+ are defined. For any two neighboring circles Co, C\ of the pattern we 



use the notations as on Fig. |11|, including the notation tp for the angle zqz%z 2 . It is not difficult to 
express the angle rp in terms of the radii ro, t\ of the circles Cq, C\ and their intersection angle <p: 



exp(iip) 



r + r\ exp(— i 
r + n exp(^ 



(7.15) 



For any circle Cq of the pattern, denote by C k , k = 1, 2, . . . , m its m neighboring circles. Let r k be 
the radius of the circle C k , denote by tp k , k = 1, 2, . . . ,m the angles ip corresponding to all pairs 
(Co,C k ), and denote by (p k , k = 1, 2, . . . , m the intersection angles of Cq with C k . Obviously, for 
geometrical reasons we have: 



\\ exp^A 



1, 



fc=i 
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therefore 

r + r k exp(-i(f) k 



TT •o^'^y-^) = i (7 . 16) 

xjl rQ + rk eX p(«0 fc ) 



Now from Lemmas n, 14 there follows 



Proposition 27 Lei G be a cellular decomposition of a topological disc (open or closed), and con- 
sider a circle pattern with the combinatorics of Q. Suppose that the intersection angles of the circles 
satisfy the condition from Proposition |7^, namely: 

• For each circle of the pattern the sum of its intersection angles 4> k , k = 1,2, ... ,m with all 
neighboring circles of the pattern satisfies 2 X^fcLi 0k = (mod 2tt) . 

Then there exists a labelling c : E(T>) i— ► S 1 = {6 € C : \8\ = 1} depending only on the intersection 
angles, such that the Euclidean radii of the circles r : V{Q*) i— ► R + satisfy the corresponding 
multiplicative hyperbolic Toda system on Q* : 

m 

-tt c k r - c k+1 r k _ 1 ^ 



Proof. Intersection angles are naturally assigned to the edges from E(Q). So, on Fig. 11 



(f>(vo,v 2 ). Comparing ( 7.16 ) with ( 7.17 ), we see that the problem in question is finding a labelling 
c : E(T>) i — ^ S 1 such that, in the notations of Fig.Q, 

-exp(i0(uo,«2)) = exp(«7r - i(/>(vo,v 2 )) = —■ 

The end of the proof is exactly the same as that of Proposition |l5|. I 

Finally, we discuss the duality for the hyperbolic shifted ratio system. As a basis for the duality 
there serves the formula 

sinh(z 2 —z\ + a 2 ) sinh(zi - z + ai) _ sinh(z 2 - z 3 + ai) sinh(z 3 - z + a 2 ) ^ ^ 

sinh(z 2 — zi — a 2 ) sinh(zi — Zo — «i) sinh(z 2 — Z3 — ai) sinh(z3 — zq — a 2 ) ' 

which is a simple consequence of (fO|), ( |7.2[ ). As usual, the latter equation allows us to introduce 
the dual system: 

Proposition 28 For any solution z : V(T>) \— > C of the hyperbolic shifted ratio system on V, the 
formula 

1 sinh (z(yx) — z(v 2 ) + a(e) ) 

Z(«2) - Z(ui) = - log ) L (7.19) 

sinh ( z(v 1) — z(v 2 ) — a(e) J 

(/or any e = (vi,v 2 ) € E(V)) correctly defines a function Z : V(V) 1— > C which is another solution 
of the same hyperbolic shifted ratio system. 



In the notation -X"(t>) = exp(2Z(v)), the latter formula reads: 

X{v 2 ) _ c(e)x(vi) - x(v 2 ) 
X(vi) x(vi) - c{t)x(v 2 )' 



(7.20) 
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In particular, on each face of T> we have: 

Xi _ CiXq- Xi X 2 _ C 2 Xi 

X x - C X X\ ' Xi 



X2 



X 2 



c\x 2 - x 3 



X 3 



c 2 x 3 - x 



X - C\X\ Al X\ - C 2 X 2 A 2 X 2 - C\X 3 A 3 X 3 

The first and the fourth of these equations together with formula (|7.9| ) yield: 

X 1 _ cix - c 2 x 2 
X~ 3 ~ 



C2X0 



(7.21) 



C 2 X — C\X 2 

This shows that, like for the additive rational Toda systems, a solution of the multiplicative hyper- 
bolic Toda system on Q determines the dual solution of the multiplicative hyperbolic Toda system 
on Q* almost uniquely (up to a constant factor, in terms of the variables x, X, or up to an additive 
constant, in terms of the variables z, Z). 

One can also consider the function w : V(T>) 1— » C which coincides with x on V{Q) and with X 
on V(Q*). Then formula (7.21) defines a system of equations on the diamond, reading on the face 
(v ,vi,v 2 ,v 3 ) as 



Wi 

w 3 



CiWq - C 2 W 2 



(7.22) 

C 2 Wq — C\W 2 

This is a well-defined equation on the quad-graph, which in the simplest possible situation (square 
lattice, constant labelling) was introduced in [H2] and therefore will be called the Hirota equation on 
T>. It appeared also in ]NQC| under the name of "discrete MKdV" . Actually, when considered on a 
single face, this equation represents nothing but the famous result of Bianchi on the permutability of 
Backlund transformations for the sine-Gordon equation. See a discussion of the role of this equation 
(on a square lattice) for the geometrical and analytic discretization of surfaces of a constant negative 
Gaussian curvature in |BP2]. 

The equation ( [7.22 ) is in the three-leg form, and restrictions of its solutions to V(Q) and to 
V(G*) are solutions of the multiplicative hyperbolic Toda systems on the corresponding graphs. 
Also this equation possesses the property of the three-dimensional consistency, which yields inte- 
grability: assign the label \i to all "vertical" edges (v,v') of D, and consider the Hirota equation 
on the "vertical" face (vq, vi, v[, v' ): 



Wn 



cw. 



CWq 



Solving for w'± in terms of w' , we find: 



CWqWq — flWoWl 



Hw'q 



CWq 
/i 



-/j-w w 1 

— CW\ 



[w' ] = L(wi,w ,c,fi)[w' ] 



where (after normalizing the determinant) 



L{wx,wq,c,h) = (fi 2 - c 2 ) 1/2 



V 



WlJ 



1/2 



-fi(w wi) 1/2 

1/2 



-1/2 _ r ?1 



Wi 

w 



(7.23) 



This is the well-known transition matrix for the Hirota equation [H2], [FV]. 

We remark also that there holds the diagram like that of Fig.^T with the replacement of the 
cross-ratio system by the hyperbolic shifted ratio system, of the additive rational Toda systems 
by the multiplicative hyperbolic Toda systems, and of the discrete KdV equation by the Hirota 
equation. 
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8 Conclusions 



The present paper is devoted to the study of discrete integrable systems on arbitrary graphs. 
We considered several fundamental examples supporting our viewpoint that integrable systems on 
quad-graphs play here a prominent role. In particular, the recently introduced systems of the Toda 
type or arbitrary graphs find their natural and simple explanation on this way. Let us indicate 
several points that were left open in this paper, and some directions for the further progress. 

1) We restricted ourselves here only to local problems of the theory of integrable systems. There- 
fore, we considered only systems on graphs having the combinatorics of a cellular decomposition 
of a topological disc. We intend to address global problems for integrable systems on graphs em- 
bedded in arbitrary Riemann surfaces in our future work. This has important applications to the 
discrete differential geometry and the discrete complex analysis. 

2) Although we have unveiled several mysteries of integrable systems here (in particular, let 
us mention the derivation of zero curvature representation from the three-dimensional consistency 
of discrete systems), there remain several phenomena which still wait for their explanation. For 
example, a deeper reasons for the existense of the three-leg form for all our equations are still 
unclear. Also the role of the duality and its interrelations with integrability calls for a better 
understanding. 

3) In all our examples the basic equation on a face of the quad-graph may be written as (cf. 
Fig. |) 

$(z iZi,Z2,Z3,a 1 ,a 2 ) = 0, 

where the function $ is a polynomial (of degree < 2), linear in each of the arguments z k . As a 
result, this equation may be solved uniquely for each of z k , the result being a linear-fractional 
(Mobius) transformation with respect to any of the remaining variables. As we have seen, the 
three-dimensional permutability of Mobius transformations of this kind lies in the heart of the 
integrability. A further important class of integrable four-point equations, with functions 4? being 



polynomials of degree < 4, linear in each of the arguments, was found in [Al] and identified as a su- 
perposition formula for Backhand transformations of the Krichever-Novikov equation. Thus, these 
integrable systems are related to an elliptic curve, where the spectral parameter of the Krichever- 



Novikov equation naturally lives |KN1 |. In a subsequent publication [AS] it will be shown that 



the corresponding equations on quad-graphs still possess the three-leg form (after parametrizing 
variables by elliptic cr-functions, Zk = a(uk))' 

a(ui - u + ai)a{u\ + u - ai) a(u 3 - u - a^VC^ + uo + «2) 
a(ui - u - ai)cr(ui + u + a\) (j(u 3 - u + a 2 )cx(u3 + u - a 2 ) 
a(u2 - uo + ct\ — O2)o-(u2 + up - a\ + a 2 ) 
a(u 2 - Uq - ol\ + a 2 )o"(^2 + uo + «i - "2) ' 

As one of the consequences, we will find the elliptic Toda system on an arbitrary cellular decom- 
position Q. In the notations of Fig.^ it reads: 

A 0~{u2k -UQ + Ctk- Clk+l)0-(u2k + Up - Ofc + Qfc+l) _ 

a(u 2 k -u -a k + a k+ i)a(u 2 k +u + a k - a k+x ) 

We remark that in the degenerate case when u-functions are replaced by sinh, the latter equation 
is again related to circle patterns with the combinatorics of Q : the radii of the circles calculated in 
the non-Euclidean (hyperbolic) geometry always deliver a solution to such a system on Q* . 
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As the simplest particular case of this system (regular square lattice, constant labelling in each 
of the two directions), we will find |AS|: 



aixj — Xj + ha) a(xj + Xj — ha) a(xj — Xj — ha) a(xj + Xj — ha) 
a(xj — Xj — ha) a(xj + Xj + ha) a(xj — Xj + ha) a(xj + Xj + ha) 

a(xj — Xj-i + ha) a(xj + Xj-\ + ha) a(xj + i — Xj — ha) a(xj + \ + Xj + ha) 
a(xj — Xj-i — ha) a(xj + Xj-i — ha) o~(xj+i — Xj + ha) a(xj + \ + Xj — ha) 

Here the index j corresponds to one lattice direction of the regular square lattice, and for the second 
lattice direction we use an index-free notation, so that the variables Xj correspond to one and the 
same value of the second lattice coordinate t € hZ, the tilded variables Xj correspond to the value 
t + h, and subtilded variables Xj correspond to the value t — h. In this interpretation, the above 
equation serves as a time discretization of the so called elliptic Toda lattice IK1JJ, 



Xj = (±j — a 2 ) 



C(xj - Xj-i) + ((xj + xj-i) - C{x j+1 - Xj) + ({x j+1 + Xj) - 2((2xj] 



Similarly, if the cellular decomposition Q is generated by the regular triangular lattice, the resulting 
discrete system may be considered as a time discretization of a novel integrable lattice which is a 



sort of a "relativistic" generalization of the elliptic Toda lattice, cf. [A2]. So, we will see that the 
root of integrability of all these discrete and continuous systems still lies in the three-dimensional 
permutability of certain Mobius transformations. 

4) In the light of above, there appears a meaningful problem of finding and classifying all 
equations $(zq, zi, Z2, £3, ai, CK2) = with analogous properties, i.e. uniquely solvable with respect 
to each argument, and leading to transformations with the three-dimensional compatibility. We 
plan to address this problem in our future research. 

5) Finally, a very intriguing problem is a generalization of the results described above to the 
case of spectral parameter belonging to higher genus algebraic curves, e.g. finding the discrete 



counterparts of the Hitchin systems in Krichever's formulation [K2]. 
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